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Dynamical systems whose symplectic structure degener- 
ates, becoming noninvertible at some points along the orbits 
are analyzed. It is shown that for systems with a finite num- 
ber of degrees of freedom, like in classical mechanics, the de- 
generacy occurs on domain walls that divide phase space into 
nonoverlapping regions each one describing a nondegenerate 
system, causally disconnected from each other. These surfaces 
are characterized by the sign of the Liouville's flux density on 
them, behaving as sources or sinks of orbits. In this latter 
case, once the system reaches the domain wall, it acquires a 
new gauge invariance and one degree of freedom is dynam- 
ically frozen, while the remaining degrees of freedom evolve 
regularly thereafter. 



I. INTRODUCTION 

A number of dynamical systems of physical interest 
possess field-dependent symplectic forms which degener- 
ate, becoming noninvertible for some particular config- 
urations. Systems as diverse as vortex interactions in 
fluids jjj, and gravitation theories in dimensions d > 4 
containing higher powers of curvature in the Lagrangian 
exhibit this feature (see e.g. Q). Models of this kind nat- 
urally arise in different contexts of current high energy 
physics, ranging from cosmology and brane worlds (g|Jj] 
to strings an M-theory 

The problem is how to describe the evolution of the 
system near a degenerate configuration and, if it could 
reach such state, how it would evolve afterwards. The 
standard hypotheses in the treatment of dynamical sys- 
tems, however, exclude the possibility that the symplec- 
tic form may have nonconstant rank throughout phase 
space, even in classical mechanics (see, e.g., |^||). 

As a first step towards understanding the general prob- 
lem, here we analyze degenerate dynamical systems in 
classical mechanics. We show that it is possible to fully 
characterize the evolution of these systems. 

It should be emphasized that this degeneracy is in- 
dependent of Poincare's classification of singularities. 
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A Poincare singularity occurs at critical points of the 
Hamiltonian, which are generically isolated, whereas the 
symplectic form degenerates on surfaces which are gener- 
ically domain walls. This kind of surfaces cannot be un- 
derstood as dense sets of Poincare singularities. Roughly 
speaking, a symplectic degeneracy is the counterpart of 
a Poincare singularity in that, in the latter the gradient 
of the Hamiltonian vanishes, whereas the former can be 
interpreted as an infinite gradient. 

The previous point can be made more explicit, by con- 
sidering the simplest example of a degenerate system, 
whose phase flow satisfies 
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with E\E 2 7^ 0, which degenerates at x 2 = 0. An equiv- 
alent formulation in the x 2 ^ region is 
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which can be viewed as a phase flow where the gradient 
of the Hamiltonian diverges as x 2 — > 0. The required 
symplectomorphism (canonical transformation) to obtain 
Eq. (^) from Eq. ([!]) is noninvertible throughout phase 
space, however. 



II. FIRST-ORDER LAGRANGIANS AND THEIR 
SYMPLECTIC FORMS 

Let us consider a system whose action is a one-form 
A, integrated over a (0 + l)-dimensional worldline em- 
bedded in a (2n + l)-dimensional spacetime of signature 
(-,+,...+), 



S[z;l,2] 



(3) 
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The field A^ is a prescribed set of 2n + 1 functions of 
the embedding coordinates which are the dynamical 
variables |10|. This action is manifestly invariant under 
reparametrizations of the worldline r — > t'(t), and dif- 
feomorphisms — > z'^(z) pd]] . Identifying the affine pa- 
rameter with the timelike embedding coordinate z° := t, 
so that z l = z l (t), the action reads 
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%;1,2]= / [A i z i +A ]dt. (4) 

The equations of motion are given by 

F ij z j + E i = 0, (5) 

where we have defined Ei = diAo — d^Ai and F^ = 
diAj — djAi. In the following, we assume Aj and Aq to 
be time- independent. 

These dynamical systems are naturally classified ac- 
cording to the rank p of the symplecticQ form Fij . Thus, 
three cases are distinguished: (A) Regular Hamiltonian 
systems, for which the symplectic form has constant max- 
imal rank, p{Fij) = 2n throughout phase space T [flo|| . 
(B) Singular or constrained Hamiltonian systems, which 
have a constant nonmaximal rank, p{Fij) — 2m < 2n 
throughout T And, (C) Degenerate systems, which 
have nonconstant rank p(Fij) throughout T. 



III. DEGENERATE SYSTEMS 

We will focus our discussion in the degenerate case (C), 
which has been traditionally left aside in the literature. 
We will assume that the zero-measure subset of T given 
by 

S = {z e T/F = 0} , (6) 

where F := det(Fy), is not dense. Thus, outside E, 
the symplectic form Fij has a constant rank 2n, and the 
dynamical structure there is described through cases (A) 
above 0. 

Under these conditions, nothing prevents the system, 
starting from a generic state for which F ^ 0, from reach- 
ing a point on E after some finite time. Having this sce- 
nario in mind, we address the following points: 

• The description of the locus of E. 

• Classification of the phase flow near E. 

• Whether E can be reached and, in that case, the fate 
of the system thereafter. 

A. Degeneracy Surfaces E 

As is well known, a skew-symmetric 2n x 2n ma- 
trix Fij(z) can be brought into the block-diagonal form 



1 The standard notion of symplectic form is usually assigned 
only to nonsingular closed two- forms. However we extend the 
term "symplectic form" to the cases (B) and (C) below. 

2 The case in which p{Fij) is less than maximal in the com- 
plement of E, is a combination of cases (B) and (C). It is 
straighforward to consider this additional complication, but 
it does not add much to deserve an extensive discussion here. 



by an orthogonal transformation. Thus the two-form 
T = jFijdz 1 A dzi can be block diagonalized in an 
open set, under a local 0(2n) coordinate transformation 

z l — > x' L (z), 

n 

T = J2fr(z)dx 2r - 1 Adx 2r . (7) 

r=l 

However, in open sets containing points of the degen- 
eracy surfaces, the Darboux-like coordinates x l cannot be 
brought into the standard canonical form, because at least 
one of the f r 's in (0) vanishes at E. Hence, further 
(finite) rescalings cannot normalize the / r 's to 1. As 
a consequence, the set E is the union of the (2n — 1)- 
dimensional surfaces 

E r = {z e r/f r (z) = 0} , 

that is, E = U" =1 E,.. 

Moreover, by virtue of the Bianchi identity {dT = 0), 
it can be shown that f T [x) depends only on the pair of 
conjugate coordinates (x 2r_1 , x 2r ). This means that the 
degeneracy surfaces are constant along the remaining co- 
ordinates. 

We assume that the / r 's are smooth Morse functions 
on the corresponding (x 2r ~ l — x 2r ) planes, which ensures 
that they possess only simple zeros except at isolated 
points; the cases where f r has zeros of higher order can 
be thought of as the merging of simple zeros. Hence, the 
level curves f r (x 2r ~ 1 , x 2r ) — divide the (x 2r ~ 1 — x 2r )~ 
plane into nonoverlapping sets and therefore, 

Lemma 1: The locus of the degeneracy surfaces E 
corresponds to a collection of domain walls, splitting the 
phase space Tinto a number of nonoverlapping regions. 

B. Characterization of the Phase Flow near E 

Generically, at a surface E r the rank p(Fij) is lowered 
by 2, and at points where k of these surfaces intersect, p 
is lowered by 2k. In a sufficiently small neighborhood of 
the surface E r , the behavior of the system is dominated 
by the dynamical variables x a — (x 2r ~ 1 — x 2r ), whose 
corresponding equations of motion can be read from Eq. 
(|)as 

€ a pf(x)± 13 = -E a , (8) 

where for simplicity, we have set r = 1, so that a and 
[3 = 1,2 and / := f\. Near a degeneracy surface E r , the 
remaining dynamical variables z a , (a — 3, 2n), behave 
like the phase space coordinates of a regular system. 

Here it is assumed that E a , remains finite and does not 
vanish on Ei (i.e., Poincare singularities are assumed to 
be located outside E), therefore, Eq. (||) implies that 
the velocity becomes tangent to the (x 1 — x 2 ) plane, be- 
cause the components x a become unbounded as the orbit 
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approaches Si, while the other components (z a ) remain 
finite. 

Due to the fact that / has a simple zero at £1, x a 
reverses its sign across the degeneracy surface. Conse- 
quently, the phase flow evolves in opposite directions on 
each side of E. Thus, in a local neighborhood of E, one 
of the following three situations occur: (a) Orbits flow 
towards E and end there, (b) the orbits originate at the 
degeneracy surface and flow away from it, or (c) the or- 
bits run parallel to E, but in opposite directions on each 
side. 

Hence, the surfaces act as sinks or sources for the orbits 
in cases (a) and (b) respectively, which naturally suggests 
a classification of the local nature of E into £( 
and YS Q > for the cases (a), (b) and (c), respectively (see 
Fig. 1). 

In all three cases there is no flux across the degeneracy 
surface, and therefore, 

Lemma 2: The regions on either side of S are causally 
disconnected and dynamically independent from each 
other. 

An immediate consequence of this, is the violation of 
Liouville's theorem at the surfaces of degeneracy. In fact, 
outside the degeneracy surfaces, the Liouville current 

f = VFz* , (9) 

is divergence-free {dij 1 = 0) by virtue of the equa- 
tions of motion and the identity di{\fFF lJ Ej) = 0, 
with F i: >Fjk — S l k . This means that Liouville's theo- 
rem holds outside E, where the dynamical behavior is 
regular. Moreover, j l has a finite limit as the system ap- 
proaches a degeneracy surface, whose only nonvanishing 
components on each side of E are 

j Q = \f\i a = sgnif^Ep . (10) 

The local character of the degeneracy surfaces E, can 
be inferred from the flux of f across a pill box enclosing 
a portion of E. The flux density $ = frii across the 
lids of the pill box is given by the projection of j l along 
the normal to the surface rij = c^i 7,1 / 2 , whose only non 
vanishing components are n a — d a \f\, that is, 

$ = -F 1 t 2 F lj E j d l F 1 ' 2 = d a fc af) E . (11) 

Note that $ is not only finite, but continuous on E. 
Therefore, 

Lemma 3: The local character of the degeneracy sur- 
faces is given by E^- 1 with 77 = s<?n($). Furthermore, 
in general, E is globally piecewise attractive (E^- 1 ) or 
repulsive (E( + )), and is of type E^ ) at the intersections 
with the surfaces 11 = {z e T/$(z) = 0} (see Fig. l.d). 
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FIG. 1. (a), (b) and (c) show the qualitative local flow in 
the neighborhood of S' - ' and E'- -', respectively. The 

global structure of degeneracy surfaces is shown in (d). 

Hence, E^ ) generically corresponds to the boundaries 
between £<") and (that is, £<W = dE^') which is 
a subset of codimension 2 in phase space. 

In the particular case, when both surfaces E and n co- 
incide on an open set, E is globally of type E^ ^. This oc- 
curs for example, if -Ei| S (o) = di(h(z*)F 1 / 2 ), whose only 
nonvanishing components are of the form E a = h(z a )d a f 
for some functions h and h 7^ p2[ . 

C. Evolution towards E ( ~^ 

The degeneracy surfaces E^ + ^and E^ - ) represent sets 
of initial and final states of the system, respectively. Con- 
figurations at a surface E^" 1 ") are unstable against small 
perturbations, and it seems unlikely that a system could 
be prepared there. On the other hand, if one considers 
the system at £(~), a small perturbation to move it away 
from the surface would require an infinite acceleration. In 
this sense, the surfaces E*- - ' represent stable final states 
for the evolution of the system, and any initial configu- 
ration sufficiently near the degeneracy surface is doomed 
to fall on it. Then, the question whether the system can 
be consistently defined on E^ - - 1 naturally arises. 

For simplicity, let us consider a system possessing a sin- 
gle surface of degeneracy which is globally of type £' A 
We will now show that when the system reaches E'~), 
two coordinates become non dynamical; the system ac- 
quires a new gauge symmetry on the degeneracy surface 
which corresponds to displacements along E(~\ and one 
degree of freedom is lost. 
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Following Dirac's approach for constrained systems 
|^3| , the action (|J) possess 2n primary constraints com- 
ing from the definition of the canonical momenta pi = 
dL/dz\ 

<f) i (z,p)=p i -A i {z,t)ttQ, (12) 

whose Poisson brackets are {cf>i, cf>j} — Fij. Outside E(~\ 
the invertibility of Fij implies that the constraints 4n 
are second class. However, at the degeneracy surface, 
the rank of Fij is reduced by two, thus, two of the (f>'s 
have vanishing Poisson brackets with the whole set of 
constraints. 

Although the constraint structure changes abruptly at 
after the system reaches this surface, its evolution 
can be described by a standard constrained system, as 
can be seen through a suitable change of basis for the 
constraints 4>i. 

Linear combinations of the form (p( a ) = e\ a j4>i, become 
first class provided e 1 ,^ are null vectors of F^. This can 
only happen at the degeneracy surface, where there are 
two of such vectors. They can be chosen so that one is 
tangent and the other is normal to the surfaces F = con- 
stant, namely, e^Fy = \djF and e l , 2 ->Fij — FijdiVF. In 
Darboux-like coordinates, their only nonvanishing com- 
ponents are e?U = e al3 dpf and e" 2 ^ = S al3 dfjf, with 
a =1,2. 

In the basis 4>i = W{a)i4>a\, with a = 3,...,2n, the 
constraint algebra reads, 

{¥>(<*),¥>(/?)} ~ ^( a )(f3)F~HdzF) 2 = fe {a){fJ) \df\ 2 , 

{^ a ,<t>b} = F ab . (13) 

^From this it is apparent that, on the surface E* the 
constraints <f( a ) have vanishing Poisson brackets, and are 
therefore candidates for first class constraints. 

In order to examine whether tp( a ) are first or second 
class at the degeneracy surface (/ = 0), it is necessary 
to compute their Poisson brackets with /. The only non 
vanishing bracket involving / is 

{/,^ (2) } = ^ 2) d Q /Hd Q /| 2 , (14) 

which cannot vanish on E because, by hypothesis, / has 
a simple zeros there. This shows that cpm is first class, 
while, (/, (/2(2)) form a conjugate pair of second class con- 
straints. 

The transformations generated by f( a ) correspond to 
Sz a = 0, and 

5 X » = {x<*,^ i p w } = ^e$ ) =?* • (15) 

Thus, the constraints <^>(i) and </?(2) generate tangent and 
normal displacements to S' — ) respectively, as expected. 
Hence, / ~ can be viewed as the gauge fixing condi- 
tion associated with the "gauge generator" ip( 2 ) ■ This is 



summarized in the following 

Lemma 4: On the degeneracy surface E^ - -*, the sys- 
tem acquires a new gauge invariance, because the second 
class constraint y>(i) becomes first class, while the num- 
ber of second class constraints (/, <£>(2),0a) remains the 
same (2n). Since each first class constraint eliminates 
one degree of freedom, we conclude that one degree of 
freedom is dynamically frozen on the degeneracy surface. 

We illustrate these results in the following examples. 

IV. EXAMPLES 

A. Simplest Degenerate System 

The simplest case of a degenerate dynamical system is 
provided by the Lagrangian 

L D = A a x a + A , (16) 

with Ai =0, Ai = xiX2, Aq = —vx\. The symplectic 
form, F a p = t a pX2, degenerates at the surface xi = , 
which is of type E^\ with r\ = sgn(v). The orbits run 
perpendicular to E^ and take a finite time to connect a 
point on the surface with a point outside. 

This example captures the essence of the behavior of 
any degenerate system in a neighborhood of a degeneracy 
surface of type E^ + ) or E'~'. In particular, the shock- 
wave solutions of Burgers' equation, 

dtu + ud x u = vd 2 u , (17) 

which is relevant in the context of turbulence, exhibit 
this behavior. These solutions are of the form 

2n 

u(x,t) = -2^J](x-z fe (t))- 1 , (18) 

k=l 

where Zk(t) are complex coordinates which come in con- 
jugate pairs and satisfy a vortex-like equation 111]. The 
corresponding equations of motion for Zk{t) can be ob- 
tained from an action of the form (Q), which for n — 1 
and z = x\ + ixi reads 

?)(£) = (»)• <19) 

whose associated Lagrangian, is precisely given by (jl^). 
This solution describes a one dimensional shock wave cen- 
tered at x = xi, with peaks at x = x\ ± x^it) of height 
=p2z/ /x-z(t), travelling outwards from x±. 

B. Coupling with a regular system 

The next example examines explicitly the fate of a de- 
generate system when it reaches a surface of type S' — 
A simple Lagrangian for which this occurs is of the form 
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L = L D (x a ) + L R (z a )-V x (x a ,z a ) 



(20) 



X 2 (t) = ±y/2ut + (X 2 (t )) 2 , 



Here. 



L D (x a ) = A a x a -H D (x a ) , 



(21) 



with a = 1, 2, is some two-dimensional degenerate system 
possessing a global surface of type £( - ) at f(x a ) = 0; 
L R (z a ) is a Regular system with Hamiltonian H R (z a ), 
and Va(^ q , z a ) is an interaction term of the form 



Va = Xf(x a )H R (z a ) 



(22) 



This coupling is chosen so that it vanishes on and 
does not change the flux density $ there, so that the 
character of the degeneracy surface does not depend on 
the coupling constant A. Note that this coupling would 
be trivial in case of nondegenerate systems. Furthermore, 
the presence of H R in the coupling implies that, besides 
the conservation of the total Hamiltonian H = Ho + 
Hji + V\ , the equations of motion 



z a = (l + Xf(x))F ab d b H, 



(23) 



give rise to a separate conservation law for H R , because 
H R = i a d a H R = 0. In turn, this implies that the re- 
maining equations of motion 



f(x)x = d a {H D + Xf(x)H R ) 



(24) 



can be integrated as an autonomous two-dimensional 
subsystem. Once these equations have been solved, and 
their solutions substituted in (^3|), it is apparent that, 
the solutions of Eqs. ( p3|) describe the same orbits as in 
the decoupled case (A = 0) but with a reparametrized 
time, 



with 



dr 
~db 



1 + Xf(x(t)) 



this 



Note that as the orbits approach the surface 
time reparametrization remains finite. 

Once the system reaches the degeneracy surface 
(f(x) — > 0), both time coordinates become identical and, 
on £( - ) , all traces of the degenerate subsystem disap- 
pear, including the information about its initial condi- 
tions x a {to). 

Thus from the moment the degeneracy surface is 
reached, the system becomes a regular one, described 
by L R (z a ), and the degrees of freedom of the degenerate 
system are forever lost. 

In order to illustrate this point, consider the degenerate 
Lagrangian given by Eq. ( |lq ) with v < 0, coupled with 
a one dimensional harmonic oscillator in the form (22). 
In that case, the total energy is £ — £r(1 + Xx 2 ) + vx\, 
where £ R is the energy of the harmonic oscillator, which 
is separately conserved. Eq. is readily integrated as 



for t < ix2i 2 t ° ))2 , and x 2 {t) = afterwards. 

Hence, the harmonic oscillator coordinates Z = z 1 +iz 2 
evolve according to 

Z(t) = Z exp(rr) , 

with \ Zq\ 2 = 2£r, where the reparametrized time is given 
by 

r = t + A[2!4 + (x 2 (t Q )) 2 ] 3 / 2 , 
of 

for t < iS2Mt t an d T = t afterwards. 



V. DISCUSSION & OVERVIEW 

The degeneracy of the symplectic form opens up the 
possibility of a violation of Liouville's theorem. In fact, 
the divergence of the current f = yF z l reads 



da 1 



diiVFF^djAo - VFF t3 8id j A a 



If Aq — —H is continuous and differentiable, the second 
term in the r.h.s. vanishes identically. However, the first 
term can give rise to a non-zero contribution, responsible 
for the jump in the flow accross S. In this sense, the 
problem we address here is the counterpart of Poincare 
classical study of singularities in the phase flow. Both 
cases correspond to different classes of possible singulari- 
ties in the phase flow, and hence, the degeneracy surfaces 
cannot be understood as a dense set of Poincare's singu- 
larities. 

It is reasonable to expect that the extension of our 
analysis to field theory would lead to the possibility that 
the symplectic form degenerates for field configurations 
where some local degrees of freedom should freeze out 
and some field components become nondynamical. In 
the case of higher dimensional gravity, this means that 
as the system reaches a degeneracy surface, some dynam- 
ical components of the metric become redundant, which 
would correspond to a sort of dynamical dimensional re- 
duction mechanism. 

The quantum mechanical analysis of this kind of de- 
generate systems, shows that there is no tunneling across 
a surface of degeneracy S, but there is a nonvanishing 
propagation amplitude between states in the bulk and 
on S fig] . These results would be relevant for the quan- 
tum Hall effect [^6|, and also for strings propa gati ng in 
a background possessing a nonconstant i?-field p7| . 
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